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The onset of convection in the form of inertial waves in a rotating fluid sphere is studied 



through a perturbation analysis in extension of earlier work by Zhang (1994) Explicit 
expressions for the dependence of the Rayleigh number on the azimuthal wavenumber 
are derived and new results for the case of a nearly thermally insulating boundary are 
obtained. 



1. Introduction 

Convection in the form of slightly modified inertial waves is a well known phenomenon 
in geophysical fluid dynamics. The analysis of the onset of convection in a horizontal 
fluid layer heated from below and rotating about a vertical axis was first done by Chan- 
drasekhar more than 50 years ago. For an account of this early work we refer to his famous 
monograph (Chandrasekhar, 1961). He found that convection sets in at high rotation 
rates in the form of modified inertial waves when the Prandtl number is less than about 
0.6 depending on the boundary conditions. Another important case in which convection 
in the form of modified inertial waves occurs is the rotating fluid sphere heated from 
within and subject to a spherically symmetric gravity field. The transition from convec- 
tion in the form of columns aligned with the axis of rotation to inertial convection in the 



form of equatorially attached modes has been demonstrated by Zhang & Busse (1987) 
In a later series of papers Zhang (1993, 1994, 1995) developed an analytical theory for the 
critical parameter values for the onset of convection based on a perturbation approach. 
The buoyancy term and viscous dissipation are introduced in the equation of motion as 
small perturbations of inviscid inertial waves and the balance of the two terms is the 
used for the determination of the critical value of the Rayleigh number. In this paper we 
extend this approach to case of a spherical boundary of low thermal conductivity on the 
one hand and to an alternate method of analysis on the other hand which will allow us to 
obtain explicit expressions for the dependence of the Rayleigh number on the azimuthal 
wave number. 



2. Mathematical Formulation of the Problem 

We consider a homogeneously heated, self-gravitating fluid sphere rotating with the 
constant angular velocity f2 about an axis fixed in space. A static state thus exists with 
the temperature distribution Ts — Tq — fir^r^ /2 and the gravity field given hy g = —^r^r 
where r is the position vector with respect to the center of the sphere and r is its length 
measured in fractions of the radius ro of the sphere. In addition to the length tq, the time 
r^/v and the temperature j^ar^ are used as scales for the dimensionless description 
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of the problem where v denotes the kinematic viscosity of the fluid and k is its thermal 
diffusivity. The density is assumed to be constant except in the gravity term where its 
temperature dependence given by a = {dg/dT)/g — const, is taken into account. The 
basic equations of motion and the heat equation for the deviation O from the static 
temperature distribution are thus given by 

dtu + rk X u -Vtt = Qr + V^u, (2.1a) 

V-M = 0, (2.16) 

Rr-u + \7^e-Pdte = 0, (2.1c) 

where the Rayleigh number R, the Coriolis parameter r and the Prandtl number P are 
defined by 

R=^, r.^, P.^. (2.2) 

We have neglected the nonlinear terms u ■ Vu and u ■ VO in equations (|2.ip since we 
restrict the attention to the problem of the onset of convection in the form of small dis- 
turbances. In the limit of high r the right hand side of equation (|2.1ap can be neglected 
and the equation for inertial waves is obtained. For the description of inertial wave solu- 
tions Uq we use the general representation in terms of poloidal and toroidal components 
for the solcnoidal field uq, 

tto = V X {Vv xr) + Vw X r. (2.3) 

By multiplying the (curl)^ and the curl of the inertial wave equation by r we obtain two 
equations for v and w, 

[dtC2 - Td^]V^v - tQw = 0, (2.4a) 
[dtC2-Td^]w + TQv^0, (2.4&) 

where dt and dip denote the partial derivatives with respect to time t and with respect to 
the angle tp of a, spherical system of coordinates r, 9, ip and where the operators C2 and 
Q are defined by 

£2 = -r^W^ + dr{r^dr), (2.5a) 
Q = rcosflV^ - (£2 + r9^)(cos6'9r - r^^ sinede). (2.56) 

General solutions in explicit form for inertial waves in rotating spheres have recently 
been obtained by Zhang et al. (2001) Here only solutions of equations (|2.4p for which 
V is symmetric with respect to the equatorial plane and does not possess a zero in its 
^-dependence are of interest since only those are connected with the preferred modes for 
onset of convection (Zhang, 1994). These modes are given by 

Vq = P™(cos0) e-xjp{imip + iujTt}f{r), wq — P™^i(cos0) exp{im(^ + iujTi\g{r), 

(2.6) 



with 



/■/ \ m m+2 / \ m+1 2im{ni + 2) 

ujQ ;^^-j7^(l ± (1 + "^("1 + 2)(2m + 3)"^)^). (2.76) 

Before considering the full problem (|2.ip we have to specify the boundary conditions. 
We shall assume a stress-free boundary with either a fixed temperature (case A) or a 
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r ■ u — r ■ V(r x u)/r'^ — and 



= case A 
drQ = case B 



at r = 1. (2.8) 



Following Zhang (1994) we use a perturbation approach for solving equations (j2.ip . 

u = Uo+Ui + ..., uj = ujQ + LUi + ... (2.9) 

The perturbation Ui consists of two parts, Ui = Ui+Ub where Ui denotes the perturbation 
of the interior flow, while Ub is the Ekman boundary flow which is required since Uq 
satisfies the first of conditions (12. 8p . but not the second one. 

After the ansatz (|2.9p has been inserted into equations (|2.1ap and (j2.1bl) we obtain 
the solvability condition for the equation (|2.1al) for Ui by multiplying it with Uq and 
averaging it over the fluid sphere, 



iuji{\uo\^) - {Qr ■ u*) + {u* ■ W^uo + «&)), 



(2.10) 



where the brackets (...) indicate the average over the fluid sphere and the * indicates the 
complex conjugate. The evaluation of the second term on the right hand side of (|2.10p 
yields 

{u* ■ V^{uq + Ub)} = ((V X u*) ■ (V X Ub))} + ^j>[ul- Vub - Mo • (r • V)ub](fS, (2.11) 

since V^mq vanishes (Zhang, 1994). Since Ub is of the order r^^/^ and vanishes outside a 
boundary layer of thickness t~^/^ only the term involving a radial derivative of Ub makes 
a contribution of the order one on the right hand side of equation (|2.11[) . This term can 
easily be evaluated because of the condition r ■ Vr x (mq + Ub)/r'^ = at the surface of 
the sphere. Using expressions (|2.6|) and ()2.7a|) we thus obtain 



{ul-V^Ub} = \ f^\K. 



^d{cos 6) m{m + 1)(2to + 1) 



2(m + 1)2 



4 + (to + 2) 



(2m + 1) 



2m + 3 



where the relationship 

»i 



IP. 



,m+l 1 2 



dcosO = 



(2m + l)(wo(m + 1)(to + 2) 
(2to+1)2 



2m + 3 



IPn'rfcose 



(2.12) 



(2.13) 



has been used. 



3. Explicit Expressions in the Limit Pt <^ 1 

The equation (|2.1cp for O can most easily be solved in the limit of vanishing tPloq. In 
this limit we obtain for 8, 

e = P^{cose)e-Kp{imip + iLUTt}h{r), (3.1) 

with 

h{r) — m{m + 1)R 



(to + 5)(m + 4) — (to + l)m (m + 3)(to + 2) — (to + l)m 

(3.2) 
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where the coefficient c is given by 

1 1 



— — case A, 

_ I (m + b)(m + 4) — (m + i )m [m + '■i)(m + 2) — [m + \)m 

"~ [lOptU JJHI^ m + 2)/m ^ ^^^^^ 

(to + 5)(m + 4) — (m + 1)to (to + 3)(to + 2) — (to + l)m 

(3.3) 

Since O is real oji must vanish according to the solvability condition (|2.10p and we obtain 
for R the final result 

R± = , ^ ' + 2to + 1 

\^(2to + 3) [(to + 1)(1 ± ^{m'^ + 4to + 3)/(2to + 3) - m]2 j 

■{2m + 9)(2to + 7)(2to + 5)^(2m + 'if /h, (3.4) 

where the two possibilities of the sign originate from the two possibilities of the sign in 
the expression ()2.76l) for tJo- The coefficient h assumes the values 

{(m + l)TO(10m + 27) case A, 

(3.5) 
(m + l)(14TO2 + 59TO + 63) case B. 

Obviously the lowest value of R is reached for to = 1 and the value i?+ for convection 
waves traveling in the retrograde direction is always lower than the value R- for the 
prograde waves. Expression (|3.4p is also of interest, however, in the case of spherical fluid 
shells when the (to = l)-mode is affected most strongly by the presence of the inner 
boundary. Convection modes corresponding to higher values of to may then become 
preferred at onset since their r-dependence decays more rapidly with distance from the 
outer boundary according to relationships (|2.7p . 



4. Solution of the Heat Equation in the General Case 

For the solution of equation (|2.1cp in the general case it is convenient to use the Green's 
function method. The Green's function G(r, a) is obtained as solution of the equation 

[dydr + ( -iworP -TO(m + 1))] G{r,a) = 5{r-a), (4.1) 

which can be solved in terms of the spherical Bessel functions jm{^-f) and ymifJ-r), 

{Gi (r, a) = A-i i,„ (ur) for ^ r < a, 

6*2 (r, a) = Ajjnif^r) + By„i{nr) for a < r ^ 1, 

where 

^J■ = V-iWorP, Ai = n (yra{p.a) - jyn{iJ-a) -^"^^f^^ ) , (4.3a, &) 

^ = -^^jm{^J■a) ^.'^')^l , B = fij,n{lJ.a). (4.3c, d) 

A solution of the equation (|2.1cP can be obtained in the form 



h{r) ^ G{r, a)m{m + l)(a™ - a"'+^)a'^da = (4.4) 

"'0 

- / G2{r,a)m{m + l){a"' - a"'+^)a^da- [ Gi{r, a)m{m + l){a"' - a"'+^)a^da. 
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Figure 1. The Rayleigh number _R as a function of tP for m — 1,2, 3, 4, 6 and 8. Results based 
on explicit expressions such as (|4.6p in the case of m = 1 (solid lines) are shown in comparison 
with the results obtained with a Galerkin numerical scheme (dotted lines for retrograde mode, 
dashed lines for prograde mode), (a) Case A, fixed temperature boundary conditions, (6) Case 
B, insulating thermal boundary conditions. 
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Evaluations of these integrals for m — 1 yield the expressions 



h{r) ^ < 



2R 



r{^? + 10) - /^r^ - 



2i? 



10(/ir cos(/ir) — sin(/ir)) ^ 
r{fi^ + 10) - /iV' 



case A, 



r2(/iCos/i — sin/i) j 
2 „2„3 (a*^ - 10)(AJrcos(/ir) - sin(/^r)) ^ ^^^^ ^ 



- (2/i cos/i — (2 — /i-^) sin/i) / 



(4.5) 



Slightly more complex expressions are obtained for m > 1. Expressions (j4.5p can now be 
used to calculate R and uji on the basis of equation (|2.10p . In the case m = 1 we obtain 



R^21{uJoTPf 1 



2 - lOSO/i"-* 



(4.6) 



350/1 ^ sin yU 



[lOpt] 



9 + 52511-^-^ 



/i cos /i — sm /i 
(7^1^ - 70 



case A, 



175/i ^)sin/i 



2/i cos /i + (/i^ — 2) sin 



case B, 



where indicates the real part of the term enclosed by {}. Expressions (|4.6p have been 
plotted together with the expressions obtained for higher values of m in figures [1^ and 
[TJj for the cases A and B, respectively. We also show numerical values by broken lines 
which have been obtained through the use of a modified version of the Galerkin method 



of Ardes et al. (1997) Because the numerical computations have been done for the finite 
value lO'^ of r the results differ slightly from those of the analytical theory. Since there are 
two values of Wq for each m, two functions R{tP) have been plotted for each to. For values 
tP of the order unity or lower, expressions p.4p are well approached and the retrograde 
mode corresponding to the positive sign in (j2.7&P yields always the lower value of R. But 
it looses its preference to the prograde mode corresponding to the negative sign in (j2.76p 
as tP becomes of the order 10 or larger depending on the particular value of to. This 
transition can be understood on the basis of the increasing difference in phase between 
Q and Ur with increasing rP. While the mode with the largest absolute value of luq is 
preferred as long as O and Ur are in phase, the mode with the minimum absolute value 
of uj becomes preferred as the phase difference increases since the latter is detrimental 
to the work done by the buoyancy force. The frequency perturbation uji usually makes 
only a small contribution to lu which tends to decrease the absolute value of lu. 

For very large values of tP the Rayleigh number R increases in proportion to {tP)^ for 
fixed to. In spite of this strong increase 8 remains of the order tP on the right hand side 
of equation (|2.1ap . The perturbation approach thus continues to be valid for r — > oo as 
long as P 1 can be assumed. For any fixed low Prandtl number, however, the onset 
of convection in the form of prograde inertial modes will be replaced with increasing r 
at some point by the onset in the form of columnar convection because the latter obeys 
an approximate asymptotic relationship for R of the form (rP)"*/^ (see, for example, 



Busse (1970) I. This second transition depends on the value of P and will occur at higher 
values of r and R for lower values of P. There is little chance that inertial convection 
occurs in the Earth's core, for instance, since P is of the order 0.03 while the usual 
estimate for t is 10^^. 
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Figure 2. The Rayleigh number i? as a function of m for 2 • 10^ ^ tP < 10"' (from bottom 
to top). The hnes are equidistant with a step of A(tP) = 400. The filled circles indicate the 
preferred values of m. The empty circles correspond to the preferred value of m in the case when 
m = 1 is not included in the competition, (a) Case A, fixed temperature boundary conditions, 
(6) Case B, insulating thermal boundary conditions. 



5. Discussion 

Since the curves R{tP, m) intersect at values of the order 10^ of rP in figures [T^ 
and HJ), a different way of plotting the results has been adopted in figure [21 Here the 
preferred value of m has been indicated by a filled circle in the case of the prograde 
inertial mode. The results of figure [2^ agree well with those of figure 4 of Zhang (1994) 
even though only an approximate method had been used for the determination of the 
Rayleigh number. Zhang neglected the (m — l)-mode and thus arrived at a different 
criterion for the preferred mode. His preferred values of m are indicated by open circles 
in figure 2a. The (m = l)-mode could indeed be suppressed by the presence of an inner 
concentric spherical boundary. A rough estimate indicates that inertial convection with 
the azimuthal wavenumber m will be affected significantly when the radius r] of the inner 
boundary exceeds a value of the order (1 — m~^). Unfortunately an analytical theory of 
inertial waves in rotating spherical fluid shells does not exist and it is thus not possible 
to extend the analysis of this paper to the case when an inner boundary is present. For 
a numerical study of inertial convection in rotating spherical fluid shells and its flnite 



amplitude properties we refer to the paper of Simitev & Busse (2003) 



The two transitions between modes of different types mentioned in the preceding sec- 



tion illuminate some of the puzzling findings of Zhang & Busse (1987) and Ardes et al. (1997) 



The transition labeled I in figure 17 of Zhang & Busse (1987) can now be clearly identi- 
fied with the transition from retrograde to prograde inertial convection. The main result 
of our analysis is that this transition depends primarily on the parameter combination 
tP with only a minor dependence on the wavenumber m. The second transition from 



8 



F. H. Busse and R. Simitev 



inertial to columnar convection can not be pinned down as well because of the lack of a 
sufficiently accurate analytical theory for thermal Rossby waves in the low Prandtl num- 
ber regime. According to the numerical results of Ardes et al. (1997) (see their figures 4 
and 5) there exists a broad transition range involving perhaps several transitions where 
the onset of convection occurs in the form of multi-cellular modes. An illumination of 
this regime should be the goal of future research. 
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